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We introduce the general formulation of a renormalization method suit- 
able to study the critical properties of non-equilibrium systems with steady- 
states: the Dynamically Driven Renormalization Group. We renormalize the 
time evolution operator by computing the rescaled time transition rate be- 
tween coarse grained states. The obtained renormalization equations are 
coupled to a stationarity condition which provides the approximate non- 
equilibrium statistical weights of steady-state configurations to be used in 
the calculations. In this way we are able to write recursion relations for the 
parameters evolution under scale change, from which we can extract numerical 
values for the critical exponents. This general framework allows the system- 
atic analysis of several models showing self-organized criticality in terms of 

usual concepts of phase transitions and critical phenomena. 
PACS numbers: 64.60.Ak, 64.60.Lx, 05.40. +j 

In the last decade non-equilibrium critical phenomena have attracted a wide interest 
in statistical physics. Critical systems are characterized by the absence of a characteristic 
lengthscale, strong fluctuations and non-analyticity of the correlation functions. Examples 
of this behavior can be found in phase transitions self-organized critical (SOC) systems 
1^, fractal growth |^ and a vast class of complex systems 0. The major source of difficulties 
in the study of non-equilibrium critical phenomena [^,0 lies in the absence of a general 
criterion, like the use of the Gibbs distribution in equilibrium systems, to assign an ensemble 
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statistical measure to a particular configuration of the system. The probability distribution 
is instead a time dependent solution of a master equation, which only in some particular 
cases becomes stationary in the long time limit. 

In this letter we present the general formalism of a real space dynamical renormalization 
group (RG) scheme for systems with a non-equilibrium critical steady-state: the Dynami- 
cally Driven Renormalization Group (DDRG). The method combines the renormalization 
of the time evolution operator with a stationarity condition which allows the calculation of 
the approximate steady-state configurations probability distribution. This coupling acts at 
each coarse graining step and therefore represents a driving for the renormalization group 
equations. For SOC systems ||8HlO|, the DDRG allows us to derive in a broader framework 



previous RG schemes |[TT|-[T3| and to formulate a more systematic approach. Here we show 



the explicit application of the DDRG to the FFM P,p!0|, which we can now study in the whole 



parameters space. Possible applications of the DDRG are not restricted to SOC models: the 
method can be used to study other equilibrium or non-equilibrium critical phenomena such 
as driven diffusive systems which to our knowledge have never been approached by 

real space RG methods. 

We consider discrete lattice models on a (i-dimensional lattice. To each site i is associated 
a variable cXi, which can assume q different values (cTi = 0, 1, ■ ■ ■ , g). A complete set a = {cTj} 
of lattice variables specifies a configuration of the system. We define (a | T(/i) | cr°) as the 
transition rate from a configuration ctq to a configuration cr in a time step t as a function of 
a set of parameters /i = {fii}. The time dependent probability distribution P{cr,t) for the 
configurations of the system, obeys the following master equation (ME) 

P{a,to + t) = J2{a iTifi) \ a^)P{a^to). (1) 

The explicit solution of the master equation is in general not available but we can extract 
the critical properties of the model by a renormalization group analysis. We coarse grain 
the system by rescaling lengths and time according to the transformation x ^ hx and 
t b^t. The renormalization transformation is constructed through an operator 7^(5", a) 



that introduces a set of coarse grained variables 5* = {Si} and rescales the lengths of the 
system [jl^. In general, 7?. is a projection operator with the properties TZ{S,a) > for any 
{Si}, {cTi}, and J2{s}'^{S,a) = 1. These properties preserve the normalization condition of 
the renormalized distribution. The explicit form of the operator TZ is defined case by case 
in the various applications of the method. Usually, it corresponds to a block transformation 
in which lattice sites are grouped together in a super-site that defines the renormalized 
variables Si by means of a majority or spanning rule. 

We subdivide the time step in intervals of the unitary time scale {to = 0) obtaining the 
coarse graining of the system as follows: 

P'iS, t') = Y: a) Y: l \ a°)P(a°, 0) (2) 

where we have included the application of the operator TZ and t' = Vt. The meaning of 
(cr I T^\ij) I 0"°) has to be defined explicitly: the simplest possibility V = N where is 
an integer number, and denotes the application of the dynamical operator times. In 
general, since we are dealing with a discrete time evolution, we have to consider T''^ as a 
convolution over different paths, chosen by an appropriate condition. The detailed definition 



of the effective operator T is reported in Ref . [|15| . By multiplying and dividing each term 
of eq. (0) by P'(S'°,0) = X]{crO} ^("S"^, (t°)P((t°, 0) and using the properties of the operator 
IZ, we get after some algebra: 

{50} Z]{a0} /<-('^'',crU)P(cr'J,0) 

which finally identifies the renormalized dynamical operator (S" | T' | S'''). In other words 
the new dynamical operator T' is the sum over all the dynamical paths of steps that from 
a starting configuration {af} lead to a configuration {(Xj}, which renormalize respectively 
in {S^} and {5*^}. The sum is weighted by the normalized statistical distribution of each 
configuration. 

We apply this scheme to systems with a steady-state described by a stationary distri- 
bution P{cr,t oo) = W{cr). For equilibrium systems the stationary distribution has the 



Gibbs form W{a) ~ exp{—f3H{a)), where H{a) is the Hamiltonian. There is not such a 
general criterion for non-equihbrium dynamical system, therefore we have developed an ap- 
proximate method to evaluate the stationary distribution to be used in the calculation of the 
renormalized master equation. The simplest approximation considers only the incoherent 
part of the stationary distribution which does not include correlations and can therefore be 
factorized. For systems characterized by a q-state variables it has the form 

W^^{a)=Yl{p.^) (4) 
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where (p^) is the average density of sites in the K-state. In this way, we have approximated 
the probability of each configuration {cXi} as the product measure of the mean field prob- 
ability to have a state in each corresponding site. The values of the densities {(Pk)} as 
a function of the parameters p are obtained by solving appropriate mean-field equations in 
the long time limit. These equations have the form of a stationarity condition 

|{(P.)} = '5,({(p.)}) = (5) 

where the operator describes the evolution of the system as a function of the dynamical 
parameters defined above. Time independent solutions of Eq. (|^) will be referred to as 
"steady-states" , although we should keep in mind that those are only the average states of 



the ensemble ||16[. In ordinary statistical systems, Eq. @ represents the thermodynamic 
equilibrium condition. For driven dynamical systems, it describes the driving of the system 
to the non-equilibrium steady-state, by means of a balance condition. 

By inserting this approximate distribution in Eq. (^), we obtain the renormalized dy- 
namical operator 

E^.o}7^(50,aO)^,(p,„) 
where the densities are calculated at each coarse graining step from the stationary condi- 
tion (Eq.|^) with the corresponding renormalized djTiamical parameters {p}. Since in this 
framework Eq. drives the RG equations acting as a feedback on the scale transformation, 
we call it the driving condition. 



The Eq.s (^,(|^) are the basic renormahzation equation from which the desired recursion 
relations are obtained. Imposing that the renormahzed operator T' has the same functional 
form of the operator T, i.e. T'(/i) = T{fi'), we obtain the rescaled parameter set = /(yu). 
This implies that the renormalized single time distribution P'{S,t') has the same functional 
form of the original distribution P{a,t). The critical behavior of the model is obtained by 
studying the fixed points /i* = fi^jJ*)- Since we are dealing with discrete evolution operators 
T, we define the time scaling factor V as the average number of steps we apply the operator 
T in order to obtain that T'(/i) = T(yu') for the coarse grained system. In this way we obtain 
a time recursion relation t' = g{n)t, or equivalently V = g{fi), from which it is possible to 
calculate the dynamical critical exponent z = logg{fi*)/logb. In this form of the DDRG, 
we take into account only the uncorrelated part of the steady-state probability distribution. 
The results obtained are not trivial because correlations in the systems are considered in 
the dynamical renormahzation of the operator T, that given a starting configuration traces 
all the possible paths leading to the renormalized final configuration. Moreover, geometrical 
correlations are treated by the operator TZ that maps the system by means of spanning con- 
ditions or majority rules. The renormalized uncorrelated part of the stationary distribution 
is evaluated from the stationary condition with renormalized parameters, thus providing an 
effective treatment of correlations. One can then improve the results by including higher 
order contributions to the unknown stationary distribution W{a) using cluster variation 
methods Ijl^. Naturally the above scheme can also be applied to equilibrium critical phe- 
neomena, where the driving condition is represented by the equilibrium mean field equations 

The DDRG is a useful tool to study the critical properties of SOC systems. In fact, 
these systems evolve spontaneously in a scale invariant stationary state. The Forest Fire 
Model is a simple automaton which has been introduced by Bak et al. ||^ as an example of 
SOC, and has been then modified by Drossel and Schwabl [jlO|. The model is defined on a 
lattice in which each site can be empty (cTj = 0), occupied by a green tree (cij = 1) or by 
a burning tree (cTj = 2). At each time step the lattice is updated as follows: i) a burning 
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tree becomes an empty site; ii) a green tree becomes a burning tree if at least one of its 
neighbors is burning; iii) a tree can grow in an empty site with probabihty p; iv) a tree 
without burning nearest neighbors becomes a burning tree with probabihty /. The model 
was first studied in the case / = for the limit of very slow tree growth {p 0). In this 
limit the critical behavior is trivial: the model shows spiral-shaped fire fronts separated by 
a diverging length ^ ~ P^'^'', where z/p ~ 1 |jl8l. In the case / > 0, the system is supposed 



to exhibit SOC under the hypothesis of a double separation of time scales: trees grow fast 
compared with the occurrence of lightnings and forest clusters burn down much faster than 
trees grow. This request is expressed by the double limit 6 = f/p and p —>■ 0. The 
critical state is characterized by a power law distribution P{s) = s^^ of the forest clusters of 
s sites (avalanches in the SOC terminology) and the average cluster radius (the correlation 
length) scales as i? ~ 6*^^^^. 

With the DDRG framework we are able to generalize a previous RG scheme 
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order to include the proper treatment of the time scaling and to study the limit / = 
(deterministic FFM). The dynamical rules of the FFM are local and the set of dynamical 
parameters, defined by /i = {/, p}, is obtained explicitly in terms of the dynamical operators 
acting on a single site, i.e. {1 \ T \ 0) = p and (2 | T | 1) = /. The relevant dynamical scales 
is defined by the burning process which occurs with probability one. We define a cell-to-site 
transformation with scale factor 6 = 2 or larger. The rules defining the cell renormalization 



operator TZ are standard geometrical spanning conditions |T9[, and their explicit form can 
be found in Ref. ||15[. The above scheme defines a finite lattice truncation on four (two) 
sites cells in d = 2 {d = 1), and denoting by an index a each cell configuration, we have 
that J2{ai} ~^ J2a- The renormalization equations that define the renormalized parameters 
can be conveniently written as 

where | a) and | a') are the cell states which renormalize respectively in | Sf) and | Si). We 
keep the subscript i since the states refer now to a single coarse grained site and not to a 



configuration of tlie system. Witli Wa we denote tlie stationary statistical weiglit of each a 
configuration. Tliis distribution is approximate following the DDRG scheme in the lowest 
order (Eq.4), in which the average steady-state densities (p^) are obtained as a function of 
/i = {f,p} from the stationary solution of dynamical mean field equations [pO|] . 

We focus our analysis in the critical region denoted by the condition / ^ p ^ 1, namely 
where the system shows critical behavior. The time scaling factor is obtained by imposing 
that the renormalized burning process occurs with probability one ((0 | T' | 2) = 1). In 
d = 1 this condition is fulfilled up to second order in / and p and gives z = 1, recovering 



the exact result of Ref. This result is due to the fact that in d = 1 there is only a 

possible way to span the cell, and consequently no proliferations are generated. In d = 2 
one has to consider the average over different paths, and new dynamical interactions are 
generated at each RG step. This is a signature that we need an approximation which 
truncates the parameter space after each iteration so that it remains closed. This is done by 
considering just the leading order in / and p in the renormalization equations, and ignoring 
any proliferations generated at each group iteration. With this scheme we obtain z = 1, 
which is not an exact result also if in good agreement with numerical simulations {z = 1.04 
2^). It is worth to remark that the DDRG allows to overcome the approximations present 



in the approach of Ref. ||12[ , where the time scaling was not properly considered because of 
the assumption of an infinite time scale separation. In addition the general scheme shown 
so far provides the inroad towards a systematic improvement of the results by introducing 
higer order correlations in the stationary distribution as discussed in Ref. |P3|| . 

Once the time scale factor is set we can write recursion relations for p and /, or equiv- 
alently 6' = x{6,p) and p' = y{6,p), evaluating the probabilities that a coarse grained cell 
grows or is struck by a lightning in steps. The driving condition and recursion relations 



derivation is long and tedious and the explicit equations are reported elsewhere |]T5[. The 
flow diagram is stable with respect to different coarse graining rules, and for d = 1 and d = 2 
we find a repulsive fixed point in 6'^ = and Pc = 0. The fixed point densities are obtained 
from the driving condition and depend on the dimensionality. In order to discuss the critical 



behavior we have to hnearize the recursion relations in the proximity of this fixed point and 
to find the relevant eigenvalues of the diagonal transformation: 

; A.= *: (8) 

6c,Pc ^ Sc,Pc 

In d = 2 the largest eigenvalue is given by Ai, which determines the leading scaling exponent 
= log 6/ log Ai = 0.7 (for 6 = 2) obtained in Ref. |jl2|. The result is in good agreement 



with numerical simulation {un = 0.6 In the limit / = the critical behavior is 

governed by the second eigenvalue A2. This eigenvalue and its relative exponent describes 
the behavior of the correlation length in the deterministic FFM. As opposed to Ai, the value 
of A2 depends on the absolute value of the time scaling factor and therefore could not 
be obtained without the DDRG formalism. The numerical value we obtain in d = 1, 2 is 
Up = log 2/ log A2 = 1.0, which is in excellent agreement with the simulation results z/p ^ 1 

Our characterization of the fiow diagram clarifies the critical nature of the model. The 
FFM is critical only for 6c = 0,pc = 0. This implies that 6,p are the control parameters 
of the model, and the critical state is reached only by a fine tuning of these parameters. 
Similar results are obtained by applying the DDRG to the sandpile model [^. These 
results allows us to clarify the meaning of SOC with respect to non equilibrium critical 
phenomena. In SOC literature it is often reported that the origin of scale invariance in 
nature lies in the absence of tuning parameter, like the critical temperature in Ising models. 
In the renormalization group language this would imply that no relevant parameters should 
be present. The situation is, however, more subtle. It has been recognized that a common 
characteristic of SOC systems is the presence of two time scales r^, the typical relaxation 
(activity) time, and the external driving time scale (often an external noise). In order to 
observe criticality the ratio T = Ta/Tf^ must be vanishingly small (T 0) p^ , p5| . With our 
approach we can recast the above concept in more formal terms. Our RG analysis shows 
the time scales ratio T is indeed the control parameter of SOC models. This parameter 
is the ratio between /, p and the burning time scale in the Forest Fire model or the sand 
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addition and the avalanche dissipation in sandpiles, but is always related to the ratio between 
different time scales. From a theoretical point of view the critical nature of SOC systems is 
not different from that of non-equilibrium phase transitions. The peculiarity of these systems 
is that close to the critical point the system is quite stable to changes of the dynamical time 
scales. In fact, the reduced control parameter which is defined as e = (T — 7^)/7^, in SOC 
system is T itself, being % = 0. This implies that if e ~ 0, even relevant changes of the 
control parameter (T nT and n < t~^) do not drive the system far from the critical 
region. Apparently the system would not be affected by changes of T, and in this sense 
SOC systems are not very sensitive to fine tuning of the control parameter. The meaning of 
SOC is then related to the widespread existence of phenomena ruled by very different time 
scales and not to the absence of relevant control parameters as often reported in literature. 

A.V. is indebted with J.M.J, van Leeuwen for very interesting discussions. The Center 
for Polymer Studies is supported by NSF. 
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